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> 1 Abstract 

in ; 

A modified-gravity theory is considered with a four-form field strength F, a variable gravitational 
coupling parameter G(F), and a standard matter action. This theory provides a concrete realization 
of the general vacuum variable q as the four-form amplitude F and allows for a study of its 



dynamics. The theory gives a fiat Friedmann-Robertson- Walker universe with rapid oscillations 
of the effective vacuum energy density (cosmological "constant"), whose amplitude drops to zero 
asymptotically. Extrapolating to the present age of the Universe, the order of magnitude of the 

b : 

average vacuum energy density agrees with the observed near-critical vacuum energy density of the 
present universe. It may even be that this type of oscillating vacuum energy density constitutes 
a significant part of the so-called cold dark matter in the standard Friedmann-Robertson- Walker 
framework. 
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I. INTRODUCTION 

In a previous article [1], we proposed to characterize a Lorentz-invariant quantum vac- 
uum by a nonzero conserved relativistic "charge" q. This approach allowed us to discuss 
the thermodynamics of the quantum vacuum, in particular, thermodynamic properties as 
stability and compressibility. We found that the vacuum energy density appears in two 
guises. 

The microscopic vacuum energy density is characterized by an ultraviolet energy scale, 
e(q) ~ Eyy. For defmiteness, we will take this energy scale Eyy to be close to the Planck 
energy scale E P i anck = a/^c 5 JGn ~ 1.22 x 10 19 GeV. The macroscopic vacuum energy den- 
sity is, however, determined by a particular thermodynamic quantity, e vac (g) = e — qde/dq, 
and it is this energy density that contributes to the effective gravitational field equations 
at low energies. For a self- sustained vacuum in full thermodynamic equilibrium and in the 
absence of matter, the effective (coarse-grained) vacuum energy density ? vac (?) is automati- 
cally nullified (without fine tuning) by the spontaneous adjustment of the vacuum variable 
q to its equilibrium value qo, so that ev ac (go) = 0. This implies that the effective cosmolog- 
ical constant A of a perfect quantum vacuum is strictly zero, which is consistent with the 
requirement of Lorentz invariance. 

The presence of thermal matter makes the vacuum state Lorentz noninvariant and leads 
to a readjustment of the variable q to a new equilibrium value, q' = qo + 5q, which shifts 
the effective vacuum energy density away from zero, ?v ac (go + 8q) ^ 0. The same happens 
with other types of perturbations that violate Lorentz invariance, such as the existence of 
a spacetime boundary or an interface. According to this approach, the present value of ? vac 
is nonzero but small because the universe is close to equilibrium and Lorentz-noninvariant 
perturbations of the quantum vacuum are small (compared with the ultraviolet scale which 
sets the microscopic energy density e). 

The situation is different for Lorentz-invariant perturbations of the vacuum, such as the 
formation of scalar condensates as discussed in Ref. jlj] or quark/gluon condensates derived 
from quantum chromodynamics (cf. Ref. js]]). In this case, the variable q shifts in such a 
way that it completely compensates the energy density of the perturbation and the effective 
cosmological constant is again zero in the new Lorentz-invariant equilibrium vacuum. 

The possible origin of the conserved vacuum charge q in the perfect Lorentz-invariant 
quantum vacuum was discussed in Ref. [l| in general terms. But a specific example was 
also given in terms of a four-form field strength F , 0]. Here, we use this 

explicit realization with a four-form field F to study the dynamics of the vacuum energy, 
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which describes the relaxation of the vacuum energy density ? vac (effective cosmological 
"constant") from its natural Planck-scale value at early times to a naturally small value 
at late times. In short, the present cosmological constant is small because the Universe 
happens to be old. 1 

The results of the present article show that, for the type of theory considered, the decay of 
6v ac is accompanied by rapid oscillations of the vacuum variable F and that the relaxation of 
? vac mimics the behavior of cold dark matter (CDM) in a standard Friedmann-Robertson- 
Walker (FRW) universe. This suggests that part of the inferred CDM may come from 
dynamic vacuum energy density and may also give a clue to the solution of the so-called 
coincidence problem 9J, namely, why the approximately constant vacuum energy density is 
precisely now of the same order as the time- dependent CDM energy density. 

These results are obtained by the following steps. In Sec. [Til a modified-gravity theory 
with a four-form field F is defined in terms of general functions for the microscopic energy 
density e(F) and variable gravitational coupling parameter G(F). In Sec. the dynamics 
of the corresponding de-Sitter universe without matter is discussed and, in Sec. HV] the 
dynamics of a flat FRW universe with matter, using simple Ansatze for the functions e(F) 
and G(F). In Sec. \V\ the approach to equilibrium in such a FRW universe is studied in 
detail and the above mentioned vacuum oscillations are established. In Sec. IVI| the main 
results are summarized. 



II. GRAVITY WITH F FIELD AND VARIABLE GRAVITATIONAL COUPLING 

Here, and in the following, the vacuum variable q is represented by a four-form field F. 2 
The corresponding action is given by a generalization of the action in which only a quadratic 
function of F is used (see, e.g., Refs. . Such a quadratic function gives rise 

to a gas-like vacuum l]. But a gas-like vacuum cannot exist in equilibrium without external 
pressure, as the equilibrium vacuum charge vanishes, go = 0. A self-sustained vacuum 
requires a more complicated function e(F) in the action, so that the equilibrium at zero 
external pressure occurs for q ^ 0. An example of an appropriate function e(F) will be 



1 An extensive but nonexhaustive list of references to research papers and reviews on the so-called "cosmo- 
logical constant problem(s)" can be found in Ref. A recent review on cosmic "dark energy" is given 
in Ref. Q. 

2 To clarify our notation, a four- form field has components F K \^ v (x) which can always be written as 
ZnXfii' \Z\g(%)\ F{ x )> m terms of the constant Levi-Civita symbol e K x^, the determinant of the metric 
g(x) = det g^^x), and a real scalar field F(x). Hence, we can simply write F if we speak about the 
four-form field. However, this scalar field F{x) is not fundamental, as will become clear later. 
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given in Sec. IIVBI 

The action is chosen as in Ref. [l] but with one important modification: Newton's constant 
Gn is replaced by a gravitational coupling parameter G which is taken to depend on the state 
of the vacuum and thus on the vacuum variable F. Such a G(F) dependence is natural and 
must, in principle, occur in the quantum vacuum. Moreover, a G(F) dependence allows the 
cosmological "constant" to change with time, which is otherwise prohibited by the Bianchi 
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identities and energy-momentum conservation 

Specifically, the action considered takes the following form (h = c = 1) 

( R 



S[A, g ,rt = - d'x V \ S \ + <F) + £- W j , (2.1a) 

F 2 = ——F kX [muF kX>iu , F kX ^ v = V [ K Ax^u] , (2.1b) 
Fe KXf , u ^/\i\, F KX ^ = Fe KX ^/^/\i\, (2.1c) 



where V M denotes a covariant derivative and a square bracket around spacetime indices 
complete antisymmetrization. The functional dependence on g has been kept implicit on 
the right-hand side of (j 2.1 aft , showing only the dependence on F = F(A,g) and ip. The 
field ip in (12. lap stands, in fact, for a generic low-energy matter field with a scalar Lagrange 
density, £ M (ip), which is assumed to be without F-field dependence (this assumption can be 
relaxed later by changing the low-energy constants in C M to F-dependent parameters). It is 
also assumed that a possible constant term A M in C M (ip) has been absorbed in e(F), so that, 
in the end, C M (tp) contains only ^-dependent terms. In this section, the low-energy fields 
are indicated by lower-case letters, namely, g^x) and ip(x), whereas the fields originating 
from the microscopic theory are indicated by upper-case letters, namely, A(x) and F(x) 



[later also $(x)]. Throughout, we use the conventions of Ref. 10], in particular, those for 
the Riemann tensor and the metric signature ( — h ++). 

The variation of the action (I2.1al) over the three-form gauge field A gives the generalized 
Maxwell equation, 



V ^ — {-dT + T^—dF— )) =0 



(2.2) 



and the variation over the metric g^ v gives the generalized Einstein equation, 

1 / 1 \ 1 dG~\F) 

R^v - -Rg^u) + — F — — — Rg^ v 



8nG(F) V 2 16tt dF 

+^(V M V,G- 1 (F) - 9yv UG-\F)) -e(F)g„ u + T™ = 0, (2.3) 
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where □ is the invariant d'Alembertian, the energy-momentum tensor of the matter 



■ fJ,U 



field if), and e the effective vacuum energy density 

e(F) = e(F) — F —^r > (2-4) 

whose precise form has been argued on thermodynamic grounds in Ref. 

At this point two remarks may be helpful. First, observe that the action (12. lap is not 
quite the one of Brans-Dicke theory 1(3, 12], as the argument of G(F) is not a fundamental 
scalar field but involves the inverse metric [needed to change the covariant tensor F K x^ u 
into a contravariant tensor F kX ^ u for the definition of F = y/F* according to (12.1bl) ]. This 
implicit metric dependence of G(F) explains the origin of the second term on the left-hand 
side of ( 12. 3p . Second, observe that the three- form gauge field A does not propagate physical 
degrees of freedom in flat spacetime 0, 8]. Still, A has gravitational effects, both classically 
in the modified-gravity theory with G = G(F) as discussed in the present article (see, in 
particular, Sec. IV Ep and quantum-mechanically already in the standard gravity theory with 



G = Gn (giving, for example, a nonvanishing gravitational trace anomaly [3j]). 
Using (I2.1cl) for F kX/1u , we obtain the Maxwell equation (I2.2p in the form 

The solution is simply 

de(F) R dG-\F) 
dF + 16n dF ^ 

with an integration constant fi. Hence, the constant /i is seen to emerge dynamically. In 
a thermodynamic equilibrium state, this constant becomes a genuine chemical potential 
corresponding to the conservation law obeyed by the vacuum "charge" q = F. Indeed, 
the integration constant // is, according to (12. 6p . thermodynamically conjugate to F in an 
equilibrium state with vanishing Ricci scalar R. 

Eliminating dG~ x jdF from (I2.3P by use of (12.61) . the generalized Einstein equation be- 
comes 

(R„ u - \ R g, u ) + i- ( V M V v G-\F)-g^UG-\F))-(t{F)- l iF) g, p +Tf y = , 

(2.7) 



8ttG(F) 



which will be used in the rest of this article, together with (12. 6p . 

Equations (12.61) and (12.71) can also be obtained if we use, instead of the original action, 
an effective action in terms of a Brans-Dicke- type scalar field $(x) with mass dimension 2, 
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setting — > F(x) afterwards. Specifically, this effective action is given by 



(2.8) 



The potential term in (12.81) contains, different from a conventional Brans-Dicke potential 
V($), a linear term, — for a constant /x of mass dimension 2. This linear term reflects 
the fact that our effective scalar field $ is not an arbitrary field but should be a conserved 
quantity, for which the constant parameter /i plays the role of a chemical potential that is 
thermodynamically conjugate to $. 

Indeed, if $ in (12.81) is replaced by a four-form field F given in terms of the three-form 
potential A, the resulting \xF term in the effective action does not contribute to the equations 
of motion (12.21) . because it is a total derivative, 

(fx ^W\iiF = -± [ d 4 x , (2.9) 

where the constant /j plays the role of a Lagrange multiplier related to the conservation of 
vacuum "charge" F (see also the discussion in Refs. j^, (|, where \i is compared with the 9 
parameter of quantum chromodynamics). 

Instead of the large microscopic energy density e(F) in the original action (12.1 al) . the 
potentially smaller macroscopic vacuum energy density pv = e {F) — fxF enters the effective 
action (12.81) . Precisely this macroscopic vacuum energy density gravitates and determines 
the cosmological term in the gravitational field equations (12.71) . 

Equations (12. 6ft and ( 12.7ft are universal: they do not depend on the particular origin of 
the vacuum field F. The F field can be replaced by any conserved variable q, as discussed in 
Ref. Observe that, for thermodynamics, the parameter \i is the quantity that is thermo- 
dynamically conjugate to q and that, for dynamics, \x plays the role of a Lagrange multiplier. 
The functions e(q) and G(q) can be considered to be phenomenological parameters in an 
effective low-energy theory (see also the general discussion in the Appendix of Ref. 131] ). 

Before we turn to the cosmological solutions of our particular F theory (12.11) . it may 
be useful to mention the connection with so-called f(R) models which have recently re- 
ceived considerable attention (see, e.g., Refs. QQ and references therein). The latter are 
purely phenomenological models, in which the linear function of the Ricci scalar R from the 
Einstein-Hilbert action term is replaced by a more general function f(R). This function 
f(R) can, in principle, be adjusted to fit the astronomical observations and to produce a 
viable cosmological model. Returning to our F theory, we can express F in terms of R by 
use of (12. 6p and substitute the resulting expression F(R) into (12. 7p . This gives an equation 
for the metric field, which is identical to the one of f(R) cosmology. (The latter result is not 
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altogether surprising as the metric F(R) model is known to be equivalent to a Brans-Dicke 
model without kinetic term 15( and the same holds for our effective action (I2.8P at the 
classical level.) In this way, the F theory introduced in this section (or, more generally, q 
theory as mentioned in the previous paragraph) may give a microscopic justification for the 
phenomenological f(R) models used in theoretical cosmology and may allow for a choice 
between different classes of model functions f(R) based on fundamental physics. 

III. DE-SITTER EXPANSION 

Let us, first, consider stationary solutions of the generalized Maxwell-Einstein equations 
from the effective action (12.81) . At this moment, we are primarily interested in the class 
of spatially flat, homogeneous, and isotropic universes. In this class, only the matter-free 
de-Sitter universe is stationary. 

The de-Sitter universe is characterized by a time-independent Hubble parameter H (that 
is, a genuine Hubble constant H), which allows us to regard this universe as a thermodynamic 
equilibrium system. Using 

R fJ ,u = -g^R, R=—12H 2 , (3.1) 
we get from (12. 6p and (12.71) two equations for the constants F and H: 

(e(F)-fiF) = ^G~\F), (3.2b) 

with /i considered given. 

Eliminating the chemical potential \x from the above equations, we find the following 
equation for F: 

where the functions e(F) and G~ 1 (F) are assumed to be known. 

The perfect quantum vacuum corresponds to H = and describes Minkowski spacetime. 
The corresponding equilibrium values F = F and /i = /i in the perfect quantum vacuum 
are determined from the following equations: 



e(F) F d€[F) 



= °' = ~dF~ 

F=F 



(3.4) 

F=F 
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which are obtained from (12. 6p and (12.71) by recalling that the perfect quantum vacuum is 
the equilibrium vacuum in the absence of matter and gravity fields (T^ = R — 0). 

If H is nonzero but small compared with the Planck energy scale, the if 2 term on the right- 
hand side of (I3.3p can be considered as a perturbation. Then, the correction 5F = F — F 
due to the expansion is given by 



F 8vr V ° J I V ° J dF 

\ F=F 

where x(-^o) is the vacuum compressibility introduced in Ref. 

F=F Q , 

Equally, the chemical potential is modified by the expansion (H ^ 0): 

3H 2 



(3.5) 



(3.6) 



de(F) 

fx = fi + dfx = 



87rG(F )F 

F=F 



(3.7) 



But, instead of fixing H, it is also possible to fix the integration constant fi. From 
(13. 2p . we then obtain the other parameters as functions of \i: H(fi), F(fi), and pv(/ i ) = 
e(/jt) — lj,F(/i). The cosmological constant A(/z) = pv(A i ) is zero for \i = /i , which corresponds 
to thermodynamic equilibrium in the absence of external pressure and expansion [internal — 
P vac (/xo) = — A(/i ) = 0]. From now on, the physical situation considered will be the one 
determined by having a fixed chemical potential \x. 

The de-Sitter universe is of interest because it is an equilibrium system and, therefore, 
may serve as the final state of a dynamic universe with matter included (see Sec. IV)) . 

IV. DYNAMICS OF A FLAT FRW UNIVERSE 
A. General equations 

The discussion of this section and the next is restricted to a spatially flat FRW universe, 



because of two reasons. The first reason is that 
observational cosmology (cf. Refs. 



latness is indicated by the data from 



20| and references therein) . The second 



reason is that flatness is a natural property of the quantum vacuum in an emergent gravity 
theory (cf. Ref. [l| and references therein). In addition, the matter energy-momentum 
tensor for the model universe is taken as that of a perfect fluid characterized by the energy 
density pm and isotropic pressure Pm- As mentioned in the previous section, the physics of 
the F field is considered to be specified by a fixed chemical potential \x. 
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For a spatially flat (k = 0) FRW universe 10J with expansion factor a(t), the homogenous 



matter has, in general, a time- dependent energy density pu(t) an d pressure -Pm(^)- Equally, 
the scalar field entering the four-form field-strength tensor (j 2. left is taken to be homogenous 
and time dependent, F KXfMU = F(t) \a(t)\ 3 e^u . 

With a time-dependent Hubble parameter H(t) = (da/dt)/a, we then have from the 
reduced Maxwell equation (12. 6ft : 

3 dG' 1 fdH tt2 \ de , , 

+ 2H 2 ) = — -n, (4.1) 



8vr dF \dt J dF 

and from the Einstein equation (12.71) : 

H 2 = ^ G p tot - EG ) (4.2a) 

2^- + 3H 2 = _ 87 rGPtot-2^G^-G^^, (4.2b) 
at at dt A 

with total energy density and pressure 

Ptot = Pv + Pm , -Ptot = -Pv + -Pm , (4.3) 
for the effective vacuum energy density 

p v (F) = -Py(F) = e(F) - pF . (4.4) 
With definition (14. 4p . the reduced Maxwell equation (14.11) can be written as 

where the overdot stands for differentiation with respect to cosmic time t. The above 
equations give automatically energy-conservation of matter, 

Pm + 3H (P M + p M ) = , (4.6) 

as should be the case for a standard matter field ip (recall that V^T^ = follows from 
the invariance of S M [g^ U) ip] under general coordinate transformations; cf. Appendix E of 
Ref. |ll|). 



B. Model for e(F) 

The equations of Sec. IIV Al allow us to study the development of the Universe from very 
small (near-Planckian) time scales to macroscopic time scales. Because the results do not 
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depend very much on the details of the functions e(F) and G(F), it is possible to choose 
the simplest functions for an exploratory investigation. The only requirements are that the 
vacuum is self- sustained [i.e., (13 .4p has a solution with nonzero Fq] and that the vacuum is 
stable [i.e., the vacuum compressibility (13.61) is positive, x(-^o) > 0]. 
A simple choice for the function e(F) is 

«(*0 = ^ -^ + ^4 - (4-7) 



2 X V ifl 3F 4 7 

where x > is a constant parameter (vacuum compressibility) and Fq the value of F in 
a particular equilibrium vacuum satisfying (13. 4p . The equilibrium value of the chemical 
potential fi in the perfect vacuum is then given by 

"■--sik- (4 - 8) 

The microscopic parameters Fq and x are presumably determined by the Planck energy 
scale, |F | ~ -Ep lanck and x ~ V e (-^o) ~ V-^pianck- Fro m Q4.8j) . we then see that |/x | ~ |-Pb|- 
Let us now rewrite our equations in microscopic (Planckian) units by introducing appropriate 
dimensionless variables /, y, u, k, h, and r: 

F = fF , y = f-l, (4.9a) 
Ai = -=-, G-\F) = k(f)\F \, (4.9b) 

H = h/^xW\, * = r^xW\, (4.9c) 

where the variable y has been introduced in anticipation of the calculations of Sec. |V] The 
corresponding normalized vacuum and matter energy densities are defined as follows: 

PV,M = , (4.10) 

X 

and Ansatz (14.71) gives 

= \ [~f + ^/ 4 ) - uf , (4.11) 

with u = u = -1/3 from (Oil . 

From the Maxwell equation (14. ip . the Friedmann equation (I4.2al) . and the matter con- 
servation equation (14.61) . we finally obtain a closed system of three ordinary differential 
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equations (ODEs) for the three dimensionless variables h, f, and tm' 
3 dk ( dh n , 2 \ dry 



3 ( \ dk df t 1 2 



T- \ h ^l- + kh ) = r v + r M , (4.12b) 
57T \ df dr J 

^ + 3h(l + w M ) r M = 0, (4.12c) 
with matter equation-of-state (EOS) parameter wm = Pm/pm- 

C. Model for G(F) 

Next, we need an appropriate Ansatz for the function G(F) or the dimensionless function 
g(f) = l/k(f) in microscopic units. There are several possible types of behavior for G(F), 
but we may reason as follows. 

It is possible that for F 2 Fq (i.e., in the gas-like vacuum) the role of the Planck scale 
is played by Ep(F) = le(-F)! 1 / 4 ~ |F| 1//2 . The gravitational coupling parameter would then 
be given by 

1 

G(F) 

This equation also gives the correct estimate for G(F) in the equilibrium vacuum: 1/G(Fq) ~ 
Ep lanck (F ) ~ |F |, according to the estimates given a few lines below (14.81) . Thus, a simple 
choice for the function G^ 1 (F) is 

G- 1 (F) = s\F\, k(f) = sf } (4.14) 

with / taken positive (in fact, / ~ 1 for F ~ Fq) and a single time-independent dimensionless 
parameter s also taken positive. 

Assuming (|4.14[) . the three ODEs (I4.12p become 

f dh , 9 \ dry 
° \ h % + f h ) = r v + ^M, (4.15b) 



El{F) ~ |F| , |F| < |F |. (4.13) 



^ + 3h (I + w M ) r M = 0, (4.15c) 

with ry = ry(f) given by (14. llj) and a single free parameter a = 3s/87r. This dimensionless 
parameter a is of order 1 if the physics of F field is solely determined by the Planck energy 
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scale (i.e., for F 2 ~ 1/x ~ Pq ~ -^pianck)- Anyway, the parameter a can be absorbed in h 
and r by the redefinition h — > h/ yfa and r — > r^fo. Henceforth, we set a = 1 in (I4.15p . so 
that there are no more free parameters except for the EOS parameter wu (taken to be time 
independent in the analysis of the next section). 



V. EQUILIBRIUM APPROACH IN A FLAT FRW UNIVERSE 
A. Equations at the equilibrium point p = pq 

Equations (I4.15ap - fl4.15cp allow us to study the evolution of the flat FRW universe to- 
wards a stationary state, if the initial universe was far away from equilibrium. The final 
state can be either the de-Sitter universe of Sec. [Ill] with Pm = and pv 7^ or the perfect 
quantum vacuum (Minkowski spacetime) with H = pu = Pv = and / = 1. Here, we 
consider the latter possibility where the system approaches one of the two perfect quantum 
vacuum states with / = 1, which correspond to either F = +|.Fo| or F = —\Fq\ for vacuum 
energy density (14.71) . 

Such an equilibrium vacuum state can be reached only if the chemical potential p corre- 
sponds to full equilibrium: p = p as given by (14. 8 p or u = u = — 1/3 in microscopic units 
f!4.9bp . Since p is an integration constant, there may be a physical reason for the special 
value p . Indeed, the starting nonequilibrium state could, in turn, be obtained by a large 
perturbation of an initial equilibrium vacuum. In this case, the integration constant would 
remember the original perfect equilibrium. (The evolution towards a de-Sitter universe for 
p 7^ po will be only briefly discussed in Sec. IV Dl ) 

In order to avoid having to consider quantum corrections to the Einstein equation, which 
typically appear near the time r ~ 1 (or t ~ tpianck = ^/-Epianck), we consider times r 3> 1, 
where the quantum corrections can be expected to be small. For these relatively large 
times, / is close to unity and we may focus on the deviation from equilibrium as given by 
the variable y defined in (14 .9 aft . 

Taking the time derivative of f!4.15bp for a = 1 and using (I4.15ap and (14. 15c|) . we obtain 



y - yh + 2(1 + y)h = -3 (l + w M ) r M , (5.1) 

where, from now on, the overdot stands for differentiation with respect to r. Next, eliminate 
the matter density r M from equations f!4.15bl) and (15. ip . in order to obtain a system of two 
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equations for the two variables y and h: 

y - yh + 2(1 + y)h = -3 (l + w M ) [yh + (1 + y)h 2 - r v ] . (5.2a) 

h + 2h 2 = ^ , (5.2b) 

where the last equation corresponds to f)4.15al) for a — 1. The dimensionless vacuum energy 
density (14. lip for the dimensionless equilibrium chemical potential u = u = — 1/3 is given 
by 

r v = ~y 2 + ~y 3 + ~y 4 , (5.3) 

which obviously vanishes in the equilibrium state y = 0. 

In order to simplify the analysis, we, first, consider matter with a nonzero time- 
independent EOS parameter, 

w M > , (5.4) 

so that the matter energy density from (14.15cj) can be neglected asymptotically, as will 
become clear later on. 



B. Vacuum oscillations 

Close to equilibrium, equations (j5.2aj) and (I5.2bl) can be linearized: 



y + 2h = 0, h = y. (5.5) 
The solution of these equations describes rapid oscillations near the equilibrium point: 

y = yo sincjr, h — h — — cos cur, (5.6a) 

UJ 

r v = -!/ 2 siii 2 wr, u 2 = 2 . (5.6b) 

The (dimensionless) oscillation period of y and h is given by 

To = 27r/ uj = 7rv^2 ~ 4.44 . (5.7) 

The corresponding oscillation period of the vacuum energy density ry is smaller by a factor 
2, so that numerically this period is given by r /2 w 2.22. Both oscillation periods will be 
manifest in the numerical results of Sec. IV Dl 
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C. Vacuum energy decay 



The neglected quadratic terms in equations (I5.2a[) and ( I5.2bl) provide the slow decay of 
the amplitudes in (15.61) . namely, the /-field oscillation amplitude yo(r), the Hubble term 
h (r), and the vacuum energy density averaged over fast oscillations (r v ) = y 2 {r)/A. 

The explicit behavior is found by expanding the functions y(r) and h(r) in powers of 1/r 
and keeping terms up to 1/r 2 : 

V= b -^ + C M , /> = ^ + ^, (5.8a) 

(5.8b) 
(5.8c) 

T T* 

where the equality sign has been used rather freely. Collecting the 1/r terms, we get 
homogeneous linear equations for b(r) and 1{t), which are actually the same as the linear 
ODEs (15. 5p with y replaced by b and h replaced by I. The solution of these equations is 
given by (15.61) with the same replacements: 



y 



T T* 

b c — b 

: - + —2- 

r r z 



I rh — I 

- + 5- 



.. b c - 2b 

y = - + — — 



b(r) = b sin lot , l(r) = l cos lot , uj 2 = 2 

UJ 



(5.9) 



where Iq and bo are numerical coefficients which ultimately determine the decay of h(r) and 
r v (r). 

In order to obtain these coefficients, we must collect the 1/r 2 terms. This leads to 
inhomogeneous linear equations for the functions m(r) and c(r). The consistency of these 
equations determines the coefficients Z an d b . It suffices to keep only the zeroth and first 
harmonics in the functions m(t) and c(t): 



m(r) = rrv ' sin lot , c(t) = c ( > + c ( ' cos lot . 
As a result, we obtain the following equations for m(r) and c(r): 



m — c 



2rh + c 



l - 2/q + - bl 



H — - (4/ — 1) cos cur . 

LO 



2l -hl-3(l + w M ) 6 2 -^o 



+ b LO (Z + 1) COS LOT . 



(5.10) 



(5.11a) 



(5.11b) 



From the consistency of these equations for the first harmonics of m and c, we obtain 



4Z - 1 = k + 1 



(5.12) 
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which gives Iq = 2/3. Similarly, we find from the zeroth harmonic of (15.11bj) 



( ^-1%) = 0, (5.13) 



which, for wm 7^ 0, gives b = 2y2/3 = ujIq. 

The above results for the coefficients l and b hold for the generic case wm > 0, as stated 
in (15.41) . For the special case wm = 0, inspection of (I4.15P shows that the same Ansatze for 
y{r) and h(r) can be used, but with the following coefficients: 

/ = 2/3, b = d M ujl , (5.14a) 



d M = 1 + <W,o [V 1 " ( 9 / 4 ) r Moo - lj , (5.14b) 

where a Kronecker delta has been employed in the expression for the damping factor dyi 
and the coefficient tmoo of the wm = asymptotic energy density tm ~ tmoo/t 2 has been 
assumed to be less than 4/9. 

Altogether, we have the following behavior of y{r), h(r), and ry{r) for r — > oo: 

2 \[2 

y ~ - dM — sin cur, (5.15a) 

3 r 

2 1 

d M cos cur) , (5.15b) 



3 r 

r v ~ t d M~2 sin2 ^ r > (5.15c) 



9 1V1 r 2 



with dimensionless frequency cu = \[2 and damping factor g?m given by (15.14bj) . This asymp 



totic solution has some remarkable properties (in a different context, the same oscillatory 
behavior of h has been found in Ref. [2l|; see also the discussion in the last paragraph of 
Sec. HT])- First, the solution depends rather weakly on the parameter w M of the matter EOS, 
which is confirmed by the numerical results of the next subsection. Second, the average 
value of the vacuum energy density decays as (rv) oc 1/r 2 and the average value of the Hub- 
ble parameter as ih) oc 1/r, while the average scale parameter increases as (a(r)) oc r 2 / 3 . 
Combined, the average vacuum energy density is found to behave as (rv) oc 1/a 3 , which is 
the same behavior as that of CDM in a standard FRW universe, as will be discussed further 
in Sec. IVEl 

D. Numerical results 

For ultrarelativistic matter (wm = 1/3), chemical potential fi = /io, and parameter a = 1, 
the numerical solution of the coupled ODEs fl4.15al) - fl4.15cl) is given in Figs. [1] and [2j The 
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behavior near r ~ 1 is only indicative, as significant quantum corrections to the classical 
Einstein equation can be expected (cf. Sec. IV Ap . Still, the numerical results show clearly 
that 

(i) the equilibrium vacuum is approached asymptotically (/ — > 1 for r — > oo); 

(ii) the FRW universe (averaged over time intervals larger than the Planck-scale oscillation 
period) does not have the expected behavior a oc r 1 / 2 for ultrarelativistic matter but 
rather a oc r 2 / 3 ; 

(iii) the same a oc r 2//3 behavior occurs if there is initially nonrelativistic matter, as demon- 
strated by Figs. [3] and H] for a relatively small initial energy density and by Fig. [5] for 
a relatively large initial energy density; 

(iv) for a chemical potential \x slightly different from the equilibrium value /io, the vacuum 
decay is displayed in Fig. [61 

The first three items of the above list of numerical results confirm the previous asymptotic 
analytic results of Sec. IV CI (these asymptotic results predict, in fact, oscillations between 
[0, 1] for the particular combinations shown on the bottom-row panels of Figs. [THH]), while 
the last item shows that, after an initial oscillating stage, the model universe approaches a 
de-Sitter stage (see, in particular, the middle panel of the second row of Fig. [6]) . 

E. Effective CDM-like behavior 

The main result of the previous two subsections can be summarized as follows: the 
oscillating vacuum energy density py[F(t)] and the corresponding oscillating gravitational 
coupling parameter G[F(t)] conspire to give the same Hubble expansion as pressureless 
matter (e.g., CDM) in a standard FRW universe with fixed gravitational coupling constant 
G = Gn- Recall that the standard behavior of the CDM energy density is given by Pcdm(^) oc 
a(t)~ 3 oc t -2 , which matches the average behavior found in (15.151) . 

The explanation is as follows. The average values of the rapidly oscillating vacuum energy 
density and vacuum pressure act as a source for the slowly varying gravitational field. The 
rapidly oscillating parts of h and y = / — 1 in the linearized equation (15.51) correspond to a 
dynamic system with Lagrangian density \{y) 2 — \^ 2 y 2 for a time- dependent homogenous 
field y = y(t). The F (or y) field has no explicit kinetic term in the action (12. lap , but 
derivatives of F appear in the generalized Einstein Eq. (12. 7p via terms with covariant 
derivatives of G~ 1 {F), which trace back to the Einstein-Hilbert-like term R/G(F) in (j 2.1 aft . 
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FIG. 1: Flat FRW universe with scale factor a(r), Hubble parameter h(r) = (da/dr)/a, 
ultrarelativistic- matter energy density Pm(t), dynamic vacuum energy density py(F) controlled 
by the vacuum variable F = F(t), and variable gravitational coupling parameter G = G{F). All 
variables are scaled to become dimensionless and are denoted by lower-case Latin letters, for ex- 
ample, ry = ry(f) and g = g(f). The specific choices for ry(f) and g(f) are given by (14, lip 
at chemical potential u = uq = —1/3 and by (|4.14p . respectively. The parameters of the cou- 
pled ODEs (I4.15P are chosen as (<j,wm) = (1> 1/3) and the boundary conditions at r = 1 are 
(a, h, f, tm) = (1, 1/2, 1/2, 1/20). The effective parameter tbb in the middle panel of the top row 
has been set to the value —3. 




FIG. 2: Same as Fig. [T]but over a longer time. 
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FIG. 3: Flat FRW universe with nonrelativistic-matter energy density tm, dynamic vacuum en- 
ergy density ry(f) controlled by the dimensionless vacuum variable /, and variable gravitational 
coupling parameter g = g(f). The parameters are chosen as (u, c, wm) = (— 1/3, 1, 0) and the 
boundary conditions at r = 1 are (a, h, /, ryi) = (1, 1/2, 1/2, 1/200). The effective parameter tbb 
in the middle top-row panel has been set to the value —3. The matter effects on the oscillation 
amplitudes of the bottom-row panels are small, as the damping factor dyi is close to 1, namely, 
dyi ~ 0.986 from (|5.14b|) with tmoo ~ 0.0125 [compare to the values of Fig. [5]. 




FIG. 4: Same as Fig. [3] but over a longer time. 
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FIG. 5: Same as Figs. [3] and [5] but with a larger initial density of nonrelativistic matter. 
Specifically, the parameters are (u, a, Wyi) = (—1/3, 1,0) and the boundary conditions at 
r = 1 are (a, h, f, tm) = (1, 1/2, 1/2, 1/2). Plotted on the bottom row are in the left panel: 
1/2 + (l/d M ) (3/4) (l/y/2)r(f - 1), in the middle panel: [(3/2) r (da/dr) /a + d M - l]/(2d M ), 
and in the right panel: (1/ffj^) (9/4) r 2 ry, for damping factor g?m ~ 0.590 from (|5.14bp with 
r Moo 0.290. 




FIG. 6: Same as Figs. [3] and 0] but with a small perturbation of the dimensionless chemical potential 
away from the equilibrium value, u = u$ + 5u = —1/3 — 1/40000, and evolved over an even longer 
time. 
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In a way, the effective Lagrange density \{jj) 2 — \oj 2 y 2 can be said to be induced by gravity. 
The pressure of this rapidly oscillating field y is now given by P = \{y) 2 — \uj 2 y 2 . In turn, 
this implies that the rapidly oscillating vacuum pressure is zero on average and that the 
main contribution of the oscillating vacuum energy density behaves effectively as cold dark 
matter. 3 

Observe that, while the F-field itself has an EOS parameter w = —1 corresponding 
to vacuum energy density, the net effect of the dampened F oscillations is to mimic the 
evolution of cold dark matter with w = in a standard flat FRW universe. As mentioned 
before, this effective EOS parameter w = is induced by the interaction of the F and gravity 
fields. 

An outstanding task is to establish the clustering properties of this type of oscillating 
vacuum energy density. A priori, we may expect the same properties as CDM, because the 
relevant astronomical length scales are very much larger than the ultraviolet length scales 
that determine the microscopic dynamics of the vacuum energy density. But surprises are, 
of course, not excluded. 



F. Extrapolation to large times 

In Sees. IV CI and IV D| we have established that the average vacuum energy density de- 
creases quadratically with cosmic time. This behavior follows, analytically, from (15. 15c|) 
and, numerically, from the bottom-right panels of Figs. [2J, HI and 

Extrapolating this evolution to the present age of the Universe (t now ~ 10 Gyr) and using 
|F | = s" 1 G~ l (F ) ~ 3/(87rGN) for a = 3s/8n = 1, the numerical value of the average 
vacuum energy density is given by 



{Pv{t now)) 



tl 



= f^Y ~ (4x 10- eV) 4 (!f^Y, (5.16) 

^now V ''now / \ 'now / 



for tpianck = V-^pianck ~ 5 x 10 44 s. The order of magnitude of the above estimate is in 
agreement with the observed vacuum energy density of the pre sent universe, which is close 



to the critical density of a standard FRW universe (cf. Refs. [17|,|20J and references therein) . 
If the behavior found had been (py) oc t~ n for an integer n^2, this agreement would be lost 
altogether. In other words, the dynamic behavior established in (15.1 5c|) is quite nontrivial. 



3 It is known that a rapidly oscillating homogeneous scalar field in a standard FRW universe corresponds 
to pressureless matter (cf. Sec. 5.4. f of Ref. [16]), but, in our case, matter plays only a secondary role 
compared with vacuum energy. Moreover, the oscillating scalar field gives an oscillating term in h(r) 
which is subleading (of order 1/t 2 ), whereas the oscillating term in (|5.15b| is already of order 1/r. 
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Let us expand on the previous remarks. For a standard flat FRW universe, the total 
energy density is, of course, always equal to the critical density p c = 3H 2 / (87tGn). But, here, 
the gravitational coupling parameter is variable, G = G(t), and there are rapid oscillations, 
so that, for example, (H) 2 ^ (H 2 ). This explains the following result for the case of a 
nonzero matter EOS parameter (wm > 0): 

fa mWW) = 5' (5 - 17) 

which is of order 1 but not exactly equal to 1. For nonrelativistic matter (% = 0), the 
right-hand side of (15.171) is multiplied by a further reduction factor = 1 — (9/4)r Moo , 
according to the results of Sec. IV CI 

Even though the order of magnitude of (15. or (15.171) appears to be relevant to the 
observed universe, the lit 2 behavior of (py) contradicts the current astronomical data on 
"cosmic acceleration" [9, is[ 19]. A related problem is the CDM-like expansion of the 
model universe, a oc t 2 ^ 3 , whereas big bang nucleosynthesis requires radiation-like expansion, 
a oc t 1 / 2 , at least for the relevant temperature range. Clearly, there are many other processes 
that intervene between the very early (Planckian) phase of the Universe and later phases 
such as the nucleosynthesis era and the present epoch. An example of a relevant process may 



be particle production (e.g., by parametric resonance ly, |22j), which can be expected to be 
effective because of the very rapid (but small-amplitude) oscillations. 4 A further possible 
source of modified vacuum energy behavior may be the change of EOS parameter wm = 1/3 
to wm = 0, which occurs when the expanding universe leaves the radiation-dominated epoch. 
Still, there is a possibility that these and other processes are only secondary effects and that 
the main mechanism of dark-energy dynamics at the early stage is the decay of vacuum 
energy density by oscillations. 

Another aspect of the large-time extrapolation concerns the variation of Newton's "con- 
stant." For the theory (12. ip and the particular Ansatz (I4.14p . the gravitational coupling 
parameter G(t) is found to relax to an equilibrium value in the following way: 

G-\t) ~ G-J (l + ^-f sin (JL\\ , (5.18) 



4 Energy exchange between dark matter and dark energy (for example, between cold dark matter and 
dynamic vacuum energy) may, in fact, be essential to explain the current epoch of cosmic acceleration; cf. 



Ref. [23(. The model considered in the present article does not allow for energy exchange between dark 
matter and dark energy, as (|4.6|) makes clear. However, it may be that the effects of such an interaction 
can be partially incorporated in our model as a small perturbation of the chemical potential away from 
the value (|4.8[) . The resulting behavior with exponential expansion setting in at large times (for small 
negative perturbations Su) is shown in Fig. [5] 
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with Co a constant of order unity, G^ a gravitational constant presumably very close to 
the Cavendish-type value for Newton's constant Gn, and t\jy = ^/x|F |/2 an ultraviolet 
timescale of the order of the corresponding Planckian time scale a/Gqo^/c 5 . The behavior 
(15.18!) . shown qualitatively by the / panels in Figs. []~H6j is very different from previous 
suggestions for the dynamics of G(t), including Dirac's original suggestion G oc 1/t (cf. 
Sec. 16.4 of Ref. Jl3]). For the present universe and the solar system in it, the gravitational 
coupling parameter (15.181) would have minuscule oscillations. Combined with the Planck- 
scale mass of the F degree of freedom (cf. the discussion in Sec. lVEj) . this would suggest that 
all solar-system experimental bounds are satisfied, but, again, surprises are not excluded. 5 



VI. CONCLUSION 

The considerations of the present article and its predecessor [1] by no means solve the 
cosmological constant problems, but may provide hints. Specifically, the new results are 

(i) a mechanism of vacuum-energy decay, which, starting from a "natural" Planck-scale 
value at very early times, leads to the correct order of magnitude (I5.16P for the present 
cosmological constant; 

(ii) the realization from result (I5.15P that a substantial part of the inferred CDM may 
come from an oscillating vacuum energy density; 

(iii) the important role of oscillations of the vacuum variable q (here, F), which drive the 
vacuum energy density oscillations responsible for the first two results. 

Expanding on the last point, another consequence of q oscillations is that they naturally lead 
to the creation of hot (ultrarelativistic) matter from the vacuum. This effective mechanism 
of energy exchange between vacuum and matter deserves further study. 
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NOTE ADDED IN PROOF 



Following up on the remarks in the last paragraph of Sec. [TT1 we have recently shown 27] 
that, close to equilibrium, the g-theory of the quantum vacuum gives rise to an effective 
/(-R)-model which belongs to the R + R 2 /M 2 class of models with a Planck-scale mass 
M ~ -Euv- We have also extended our analysis to a quantum vacuum containing several 
conserved g-fields, which allows for the coexistence of different vacua. 
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